Superconducting States in pseudo-Landau Levels of Strained Graphene 
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We describe the formation of superconducting states in graphene in the presence of pseudo-Landau 
levels induced by strain, when time reversal symmetry is preserved. Due to the incompressibilities 
associated with completely filled pseudo-Landau levels, superconductivity in strained graphene is 
quantum critical at integer filling factors, whereas at partial fillings superconductivity survives at 
weak coupling. In this limit, the critical temperature scales linearly with the coupling strength 
and can be orders of magnitude larger than in conventional weak coupling superconductors. We 
argue that superconductivity can be induced by electron-phonon coupling and that the transition 
temperature can be controlled with the amount of strain and with the filling fraction of the Landau 
levels. 
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Graphene is a single atomic sheet of carbon with elec- 
tronic excitations that behave as massless Dirac quasi- 
particles P, [![ ■ In general, graphene seems to be insen- 
sitive to electronic many body instabilities except in 
the quantum Hall regime where fractional quantum 
Hall states [^-Q have been observed. We claim that one 
promising route to induce intrinsic superconductivity in 
graphene is to reconstruct the electronic density of states 
(DOS) into a discrete spectrum of Landau levels (LLs) 
with the application of strong strain fields. Strain con- 
figurations in a specific engineered form Q mimic the 
application of strong uniform magnetic fields that can be 
as large as SOOT ^] , but preserve time reversal symmetry 
(TRS). Quantum Hall states induced by pseudomagnetic 
fields have been conjectured to give rise to topological or- 
der in strained graphene with spontaneously broken TRS 

In this letter, we describe the formation of intrasublat- 
tice TRS spin singlet states which occupy the LLs pro- 
duced by elastic deformations in graphene. Since the 
overall wavefunction is anti-symmetric, the spin singlet 
wavefunctions are even under valley exchange, and are 
robust against backscattering [l2| , unlike the correspond- 
ing triplet states, which break inversion symmetry fbOI. 
We show that at integer filling factors where the normal 
state becomes incompressible due to Pauli blocking, the 
superconducting order parameter has a quantum critical 
point at the mean field level, A cx |a; — xd^^^, where 
X c IS a LL dependent critical coupling, in contrast with 
the unstrained case, which is quantum critical only at 
the neutrality point jl3nl5j . Near complete filling of the 
LLs, a sequence of quantum critical points can be exper- 
imentally accessed by controlling the filling factor in the 
weak coupling regime, x ^ Xc- At partial filling of the 
LLs, we show that the zero temperature gap A oc x has a 
linear scaling with coupling and strain in the weak cou- 
pling limit a: <C Xc, and can be made orders of magnitude 
larger than in conventional weak coupling superconduc- 
tors, where A cx e~^/^ . We identify experimental signa- 



tures for this state, and we propose that in the presence 
of substrates that screen Coulomb interactions at length 
scales larger than the magnetic length, superconductivity 
can be triggered by conventional electron-phonon cou- 
pling. 

In the continuum description of the problem, the low 
energy Hamiltonian of strained graphene is 

■Ho = ^ / dx^'t „(x) [v(-iV-f aA^)-^^-^] ^...^(x) 

(1) 

where a =t, i is the spin index, a = ± indexes the two 
valleys, (T^ = {aax, (Jy) is a vector of Pauli matrices, /i is 
the chemical potential away from half filhng, v ~ 6eVA 
is the Fermi velocity, = {ipa,+,,^Pb,+ ,i>a,-,'^b,-)a is 
a 4 component spinor in the (a, b) sublattice pseudospin 
and in the two valleys and A^ is the pseudo vector po- 
tential, which couples to the electrons as a magnetic field 
pointing in opposite directions in the two different val- 
leys, preserving time reversal invariance. From now on, 
we assume strain configurations which produce uniform 
pseudo-magnetic fields, -6^,8]. 

In the Landau gauge, where A^ = (— B^y, 0), and Bg 
is the pseudomagnetic field, the electronic wavefunction 
takes the form "^k,tT{x, y) = exp(ifcx) Qcriu), where Qaiy) 
is the eigenspinor of a ID Hamiltonian. This Hamil- 
tonian can be expressed in terms of ladder operators 
of the ID harmonic oscillator, a = {£_ + d^) /\/2, = 
{£, — d^) /\/2, where ^ = £5/0 — y/^^isa dimensionless 
variable related to the valley dependent guiding center 
X = —M\, with = y^JeBs (restoring K) the effec- 
tive magnetic length, and e the electron charge. In what 
follows, we define the valley dependent operator X'(^), 



v{-iS/ + As) ■ a = V2- 

which takes the form — iV 
opposite valley, with ^ = Ib^ ^ 



^(0, (2) 



-A,) 

vIIb. 



a = -ViX) in the 
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In the presence of an effective attractive poten- 
tial U that stabiHzes the superconducting state, the 
Bogohubov-deGenncs (BdG) Hamiltonian is Hbg — 
J dx$t(x)HsG$(x), 



"Hbg = 



Ho (A, 



A 



A* -THo{A,)T- 



(3) 



where ?^(As) = ® — ® — fiaovo, is the 
normal state Hamiltonian of strained graphene written 
in valley and sublattice spaces, = {uq ± Vz)/'2. are 
projectors in the ± valley spaces, with Pauli matrices 
Vi{i = x,y,z), and $ = {'^k,-[,^ _j. ^) is the 8 compo- 
nent spinor in the Nambu space, with Pauli matrices 
Ti (i — x,y,z). The off diagonal term. A, is a pair- 
ing matrix that describes the formation of Cooper pairs 
A^ = C/tr(^'j^ ^Avf-fc^-o-). In strained graphene, the 

time reversal symmetry operation T'Ho{A.s)T^^ leaves 
the Hamiltonian invariant under an additional exchange 
between valleys, in contrast with the case of conventional 
magnetic fields, which explicitly break TRS [3]. 

In the intra-sublattice s-wave pairing state, which cor- 
responds to the pairing matrix A — AaQi>x, the eigen- 
vector problem HBG'f (s;,^) — E^{x,£^) can be solved by 
decomposing Hamiltonian ^ into two equivalent copies 
of 4 X 4 BdG Hamiltonians in pscudospin and Nambu 
spaces. In the reduced Nambu basis $ = (^^,+, 4*1 _ 



n 



BG 



ViO A 

A* + M 



(4) 



Fixing the gauge of the gap A to be real, the eigenvalue 
problem [E — Hbg)^ = is equivalent to 

M^=iE + 'H'bg){E - Ubg)^ = 

where = CHbgC^^ = (X> + /i) (g) rg + Ati is the 

charge conjugated BdG Hamiltonian (/i — > — /i). When 
the matrix M is applied in a Nambu state $ where 

^i^'' ={ip\N\-i.sip\N\)cr.a, with N the index of the Lan- 
dau levels, and s{N) = sgn(iV) accounts for the two 
branches of LLs in the conduction and valence bands, 
M can be cast in the form = {E + -H-)^, 
with 'H± ~ {sujc\/\N\ ± /i)r3 + Ari, and ujc = \/2v/£b- 
%- is equivalent to Hamiltonian and gives the energy 
spectrum 



± En = ±[(swc 



M- 




Figure 1: Dependence of the chemical potential fi in units 
of uic with the filling factor v. Red (light) curve: normal 
state A = 0. Solid black curves: superconducting state, with 
A/uic ranging from (red hne) to 1. All curves are plotted 
at T/ujc = 0.02. Inset: detail of the zero LL (n = 0). The red 
curve is analytically described by Eq. (jlOp . while solid ones 
at A/ojc <. 1 are described by Eq. @ in the T — > hmit. 



Tl/V2[l T (sWcvWl ~ tJ^)/ENY''^- 4>N{.i) denotes con- 
ventional LL wavefunctions, with = 0. In the zero 
LL, the Cooper pairs occupy only one sublattice, explic- 
itly breaking the Z2 sublattice symmetry of graphene. As 
anticipated, the BdG quasiparticle spectrum is discrete 
and can be indexed by the LL index N . The supercon- 
ducting ground state is given by 17 1 



1*0) = n iy-N + VnC^n,XA-A 



N,-X,i,+ 



|0), (7) 



(5) 



where un = /3+,t-' ""n = ^^^^ c\ x,a,a ^re 

fermionic creation operators of the relativistic LLs. This 
wavefunction describes intrasublattice pairing across op- 
posite valleys, within the same LL. 

The discontinuous behavior of the chemical potential 
with the pseudomagnetic field can be calculated by fixing 
the total number of particles in the system. Although 
the ground state wave function does not conserve the 
number of particles, the distribution is sharply peaked 
around the average in the thermodynamic limit [lg | , Af = 
gN^EN^-oo [u%fiEN) + v%f{-EN)] , where f{E) = 
(l-|-e^/"^)~^ is the Fermi distribution, g = 4 is the valley 
and spin degeneracy, and — A/{2tt£\) is the number 
of flux quanta for a total area A, which sets the LLs 
degeneracy. In the low temperature and weak coupling 
regime T, A ^ w/^s, where the deep energy states N < n 
are fully occupied, with n the highest occupied LL, the 
constraint becomes 



with eigenstates given by 



N ( (l>\N\-l{0 



Akx 



(6) 



for the states (cr, a) 

dN ^ 



S(t)\N\{0 

(t, — ) and where 



2{iy/g -n) = -[{sLj,,^\n\ - ^)/£;„] tanhK/(2T)], (8) 

where u = N/N^ — g{Nf^ + 1/2) is the filling factor, and 
Nty = (D/lJc)'^ > is an ultraviolet cutoff that regular- 
izes the number of negative energy states, where D 6eV 
is the bandwidth. In particular, at T = 0, the chemical 



potential 



/i(0,i.) 



A{iy- gn) 



(9) 

remains pinned to the n-th LL when half filled {v — gn) 
for small A, and shows a power law divergence when 
the highest occupied LL is completely filled, at integer 
fillings v = g{n ± 1/2), indicating an incompressibility. 
In the opposite regime, when T, A > w/^s, the system 
crosses over to the usual Fermi liquid behavior when the 
electrons have multiple transitions between different LL. 
In the normal state (A = 0), as expected, the chemical 
potential 



/x(T, v) — suJc\/\n\ + Tin 



gin 



(10) 



is also pinned to the energy of the highest occupied LL 
at partial fillings v — gn, and shows a logarithmic diver- 
gence at integer fillings v = g{n ± 1/2), when the corre- 
sponding LL becomes incompressible. When all LLs are 
taken into account, as shown in the red (light) curve in 
Fig. 1, at integer fillings vi{n) = g{n + i), the chemical 
potential for the normal state does not diverge but sits 
half way between the LLs in the the zero temperature 
limit, h{n) = ^[0, j/i(n)]/a;c = (v/|nTT| + v/H)/2, un- 
like the T ^ Q limit at partial filling, where ^(0, v) = 

In Fig. 1, we compute numerically the contribution 
of all LLs to the chemical potential at fixed temperature 
for different values of A. At v — the system is half 
filled and the chemical potential is always pinned to the 
energy of the zero LL. When A ~ v/is, the chemical 
potential ijl{v) detaches from the discrete profile of LL 
plateaus, and becomes a smooth function. Eq. Q and 
pop describe analytically the numerical curves shown in 
the inset of Fig. 1 when A/wc <C 1. 

At the mean field level, the free energy of the super- 
conducting state is F = ~TgN^Y,^=±H'N=-oo^'^i^ + 



-lEt,IT 



) - A|A|^/C/, where A is the total area normal- 



ized by the size of the unit cell, 
energy gives the gap equation 



Minimization of the free 



\ = -{Ul2)gN^ tan\i[EN{T,v)/{2T)]/EN, 

N=-oo 

(11) 

where = 3V3a^/ (47ri?^) is the number of flux quanta 
per unit cell, with a ~ 1.42 A the lattice spacing. Defining 
X = \U\gN^/ijJc oc \U\/tB as the dimensionless coupling 
parameter that controls the strength of interactions and 
strain, at half filling {v — Q), the zero temperature gap 
in the weak coupling regime <^ v/£b is 

A<^'\0)=\vx/iV2£B)\ /[l-a(l/2)x], (12) 
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Figure 2: Critical temperature Tc/uJc vs. the dimensionless 
coupling strength x = \U\gN^lB / {y/^v)- a) red (light) curve: 
v = 2, where the transition is quantum critical below Xq ~ 
0.025. Solid black curves: v = 0, 8, 24 and v = 40, from right 
to left, h) V = 0, 1.2, 1.6, 1.98, 1.9998 and 2, from left to 
right, c) u = 4, 5.2, 5.6 5.98, 5.9998 and 6, from left to right. 
At partial filling of the LL, Tc oc a; in the a; — >■ limit. 



where Ca(^) = X]jv=i 1/V^ is the zeta function regu- 
larized by an ultraviolet cut-off. The ratio between the 
critical temperature and the zero temperature gap at half 
filling is a universal number, 2Tc — A'^"^(0). In the weak 
coupling limit, ~ ^f^vxjlB oc Bs has a linear scaling 
with the coupling and with the amount of strain [l^ . 
This scaling contrasts with the case of conventional weak 
coupling superconductors, where Tf. oc exp(— l/x) de- 
creases exponentially with the effective coupling. As the 
coupling X becomes larger, the system eventually crosses 
over to the strong coupling regime, when ^ v/£b, as 
shown in Fig. 2. In the critical regime, when A/Tc <C 1, 
the gap at 1/ = is given by 

A(o) (T) = 2i {v/(.b)"^ v/1 - r/r„ /[C (3/2) TJ ^ (13) 

at weak coupling, and scales with A*^*^) cx bI^'^, where 
C(3/2) K, 2.61 is a zeta function. In Fig. 2a, we show the 
dependence of the critical temperature with the coupling 
X for different filling factors. The red curve is the phase 
transition for v = 2, which is quantum critical. 

At integer filhng factors viin) = 5(n + |), the 
normal state becomes incompressible at temperatures 
much smaller that the LL spacing due to Pauli block- 
ing, inhibiting charge density fluctuations required to 
preserve the electromagnetic gauge invariance in a 
superconductor [20| . In that case, the phase transition 
becomes quantum critical, and the emergence of super- 
conductivity requires a critical coupling to overcome the 
electronic incompressibility of the normal state. At those 
integer fillings, the zero temperature gap is 



A(''')(0) = 2^n{v/iB)^/x/x'i^ - 1, (14) 



where 




\N\-h{n)\-^ 



(15) 



4 




Figure 3: Phase diagram Tc/u)c as a function of the filhng 
factor u for different couphng strengths, x (soUd curves). Blue 
(dark) region: x < x\ Ki 0.022. Grey region: < a: < ajg ~ 
0.025. Light region: x > xfj. The critical temperature drops 
to zero at integer filling factors iA.{n) — g{n + ^) whenever 
X < x'fi^, where is the critical coupling of the t'i(n) state. 



is the quantum critical coupling of the = + ^) 
state, when the n-th LL is completely filled, = 
Sw=-oo IsVl-^l — is a constant and h(n) was de- 

fined below Eq. (fTU|) . For a magnetic length of £b ^ 30 A, 
which corresponds to a pseudomagmetic field Bg ^ lOOT, 
we have A^a ^ 400. At n = 0, the critical coupling is 
Xq « 0.025, as depicted by the arrow in Fig. 2b for the 
red curve, crossing over to a smooth transition for par- 
tial filling factors, as shown in the black curves in the 
same panel, where < < 2. For n = 1, the critical 
couphng drops to xl ~ 0.022 (Fig. 2c), with the sohd 
curves in the same panel indicating a crossover to partial 
filling factors in the range 4 < < 6. For higher LL, x^^ 
decreases further as n becomes large. 

For fixed coupling strength, the critical temperature 
evolves with v as a series of lobes and drops to zero at 



integer filling factors i'i{n) — g{n + ^), whenever x < x^, 
as shown in Fig. 3. The solid blue area in Fig. 3 depicts 
the region x < xl, which is quantum critical for the first 
two Landau levels, while the gray region xl < x < Xq 
is quantum critical in the zero LL only. When x grows 
larger than x"^, the system undergoes a quantum phase 
transition in the state v = vi (n) , when superconductivity 
emerges at that filling factor. For v ~ ^'i('t-), the critical 
temperature scales as 



Tc{x, ly) (x\iy - ^'i(n)|'^ 



(16) 



for X < x"^, with possible logarithmic behavior, where 
5 0.2(3) is the exponent numerically extracted for the 
v — ±2 states. This behavior may lead to the experi- 
mental observation of quantum criticality in graphene by 
controlling the filling factor of the LL in the weak cou- 
pling limit X x'^. In two dimensions, the critical tem- 
perature sets the onset of Cooper pair formation, while 
long range phase coherence is lost above the Kosterlitz- 
Thouless transition temperature Tkt < [2l|, where 
pairs of vortices and anti-vortices unbind. Those effects 



will likely be relevant for transport and will be discussed 
elsewhere. 

Besides transport measurements in long graphene 
junctions, where a supercurrent is expected to fiow along 



the edges [22|, one experimental signature of supercon- 



ductivity in strained graphene is the specific heat at fixed 
volume, Cv = ~T{d'^F/dT'^)v- At low temperature, 
T ^ Tc ^ v/£b, the specific heat of the ly = state 
is CviT) ^ 2giV0A2(O)e-^(o)/'^/T2. In the quantum 
limit, uJc/Tc S> 1, the specific heat jump at the phase 
transition normalized by the specific heat in the normal 
side for 1/ is AC/C„ = a;ce"<=/^<=/[16C(3/2)Tj which 
is non-universal. In the weak coupling regime x <^ Xq, 
where Tc xlOc/A, 



{AC/a,)(x) 



- p4/x 



/[4xC(3/2)] 



(17) 



becomes exponentially large as drops to zero. This 
feature is a signature of this state, and contrasts both 
with the specific heat jump expected for Dirac fermion 
superconductivity in unstrained graphene at half filling 
(AC/C„ ~ 0.35) [l3| and in weak^coupling superconduc- 
tors in general {l^C/Cn ~ 1.43) [isj, which are universal 
constants. 

Although strong Coulomb interactions inhibit super- 
conductivity and can give rise to incompressible states at 
fractional filling factors Q , a condensate can be induced 
by phonons in the presence of substrates that screen the 
electronic repulsion at length scales larger than £b- The 
analysis of scanning tunneling spectroscopy experiments 
[2^ in graphene on SiC for magnetic fields around 5T, 
when the LLs are well defined, indicate that the effec- 
tive momentum independent electron phonon vertex in 
graphene is go ~ O.lcV for an Einstein phonon mode 



at the typical frequency Wph ~ 0.2eV [23|, |2j|. This 
mode alone (i?2g phonon) leads to an effective attrac- 
tion U ^ — 25Q/a;ph w —0.1 eV. For a magnetic length of 
20A, which corresponds to 2 x 10~^ flux quanta per unit 
cell, a net attractive coupling of that order results in a 
dimensionless coupling x ^ 0.003 and a critical temper- 
ature Tc ~ 8K at the v = state. 

A significant enhancement of the electron-phonon cou- 
pling, and as a result Tc, can be achieved for instance by 
coating graphene with ionic crystals and alkaline metals 
such as K, which is known to form a stable crystal on top 
of graphene |25l427| . This mechanism can lead to mea- 
surable transition temperatures in the regime where the 
broadening of the highest occupied LL due to disorder 
effects is small compared to the level spacing. Our anal- 
ysis shows that the spin singlet states are robust, and 
present a sequence of quantum critical points, which can 
be experimentally acessed by tuning the filling factor of 
the LLs in the weak coupling limit of the problem. 
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